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Self-shunted aluminum Josephson junctions with high-transparency barriers were fabricated us- 
ing the shadow-evaporation technique and measured at low temperatures, T ~ 25 mK. Due to 
high junction transparency, the /^/-characteristics were found to be of only small hysteresis with 
retrapping-to-switching current ratio of up to 80 %. The observed critical currents were close to 
the Ambegaokar-Baratoff values (up to 80 — 100 %). Good barrier quality was confirmed by the low 
subgap leakage currents in the quasiparticle branches, which makes the self-shunted Al junctions 
promising for application in integrated RSFQ-qubit circuitry. 



Operation of superconducting electronic devices often 
implies an overdamped regime of Josephson dynamics 
. An important example is given by a family of Rapid 
Single Flux Quantum (RSFQ) devices (see, e.g., a re- 
view Ref. 0) traditionally based on resistively shunted 
junctions. A different type of shunting is however pre- 
ferred, when using RSFQ-circuitry as control electron- 
ics for Josephson qubits |3j. For keeping low the overall 
decoherence, low- noise behaviour of the control RSFQ- 
circuit is required at qubit frequencies of the order of 
ten GHz (see, e.g., the analysis of decoherence of the 
flux qubit 0]). Due to large Niquist- Johnson fluctua- 
tions at low frequencies, this requirement can hardly be 
met with a standard low-ohmic linear shunting. Recently, 
several realizations of frequency-dependent damping for 
Josephson junctions have been proposed [fj, [fj, @, In 
Refs. [1,0], superconductor-insulator-normal metal (SIN) 
contacts were suggested as nonlinear shunts for Joseph- 
son junctions and their operation was analyzed in simple 
SFQ networks. A detailed study of subgap conductiv- 
ity in highly-transparent SIN-junctions was presented in 
Ref. 0- Another solution on the base of RC-shunting 
was considered by Hassel et al. , to be implemented in 
form of on-chip capacitors, added in series to the stan- 
dard ohmic shunts. 

In this Letter, we suggest an alternative approach on 
the basis of Al self-shunted Josephson junctions, making 
use of their highly nonlinear quasiparticle branches. The 
power spectrum Pi qp (o;) of current fluctuations of such a 
junction in the superconducting (S) state, (V) = 0, 

P TqP H = {e/ir)I qp {hu;/e) coth(huj/2k B T), (1) 

depends on the shape of its quasiparticle branch I qp {V) 
with the current suppressed at low temperatures, 
fceT <C A, in the subgap voltage range, V < V s = 2A/e. 
For Al junctions with the superconducting energy gap 
A ~ 200 /xeV, the gap frequency, uig = 2A/H ~ 27r x 
10 n s _1 , far exceeds the typical qubit frequencies, thus 
enabling low decoherence on a time scale of qubit opera- 
tion. 

The fabrication routines for small A1/A10 X /A1 tun- 
nel junctions are well developed, basing on the shadow 
evaporation technique E3 |. In the early experiments, 



this technique has also been implemented for submicron 
self-shunted tunnel junctions composed of Pb-In/Pb [ill ]. 
Moreover, in the last years, it is widely applied for fabri- 
cation of Al-based qub its (see, e.g., the review on Joseph- 
son qubits in Ref. [l2| and references therein). The lat- 
ter fact enables, in principle, full on-chip integration of 
Al-RSFQ and Al-qubit circuitry, both operating in sub- 
Kelvin temperature range. 

The intrinsic damping properties of an autonomous 
Josephson tunnel junction crucially depend on the bar- 
rier transparency, i.e., on the specific barrier resistance 
p = R^A, where i?N is the normal junction resistance 
and A is the junction area. Using for the critical cur- 
rent its Ambegaokar-Baratoff (AB) value at T = 0, 
I c ~ I" 3 = 7rA/2ei?N, one can express the Stewart- 
McCumber damping parameter, (3 C — 2eI c Rf i C/fi — 
■nApc/Ti as a function of p, A, and the specific junc- 
tion capacitance c = C/A. The overdamping condition, 
(3 C < 1, thus imposes an upper boundary on the barrier 
transparency, p < h/(nAc). 

In experiment, non-hysteretic TV-curves [ill 0] and 
high-frequency operation of an RSFQ circuit [lj| built on 
the self-shunted Nb/A10 x /Nb junctions with ultrahigh- 
transparency barriers, p ~ 1 ft x /im 2 (corresponding to 
the high critical current density J c ~ 1 — 2mA//im 2 ), 
have been convincingly demonstrated. However the bar- 
rier quality was found to be not so high, due to noticeable 
density of pin-hole defects 0,0] and, as a consequence, 
the considerable subgap leakage currents 0, 0, 0, 0] . 

An important advantage of Al-based over Nb-based 
junctions for the specific RSFQ-Qubit applications turns 
out to be a 7 — 8 times smaller value of A « 180 /icV 
(corresponds to the zero-temperature BCS value Ao = 
1.76/cbTc for T c = 1.2 K), which allows non-hysteretic 
behaviour of more opaque junctions with, typically, p < 
30 f2 x pm 2 . This decreases the risk of pin-hole defects 
and dramatically weakens the subharmonic gap structure 
of multiple Andreev reflection |l5l ]. In the reported here 
experiments with Al junctions, we found the values of I c 
to be slightly below their AB values. We also show that 
the nonlinearity parameter r\ = [G(0)i?N] _1 ~ 10 — 30, 
where G(0) is the zero-bias conductance, substantially 
exceeds the values r\ ~ 1 measured in high-transparency 
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FIG. 1: SEM-micrograph of the two-junction SQUID. The 
junction is formed between the first and the third layers of 
Al, whereas the second layer (middle shadow of the mask) 
is not used in this structure. The /U-curves were measured 
using a four-point configuration. 



Nb-junctions [II El El. 

In order to enable on-chip comparison between junc- 
tions with different parameters, the sample was designed 
to include three consecutive Al layers ("shadows"), 30, 
40 and 20 nm thick. The layers were deposited in situ 
at three different angles, —22°, +11° and +22°, respec- 
tively, through the same PMMA/Ge/Copolymer mask, 
0.1/0.05/1.2 /on-thick, with suspended bridges. The first 
(the rightmost) shadow of Al was weakly oxidized at 
Po 2 = 0.1 Pa for 5 min to form a thin tunnel barrier. 
A set of junctions with two different transparencies was 
formed in overlaps between the first and either the sec- 
ond or, after prolonged oxidation of the first Al shadow 
by adsorbed oxygen, the third layer (cf. similar de- 
sign for Al/Cu/Cu SIN-junctions in Ref. 7]). The re- 
sults are reported for two- and four-junction SQUIDs 
with a total junction area A w 2 x (0.25 x 1) /tm 2 and 
A « 4 x (0.25 x 0.5) /mi 2 , respectively; the SQUID loop 
area was S ~ (1 x 1) /im 2 . A typical SQUID-device is 
shown in Fig. ^ We found an on- wafer spread of tun- 
nel resistances of less than 10 — 20 %, indicating good 
junction uniformity. 

The /V-characteristics (see an example in Fig. Ufa)) 
were measured at T ss 25 mK in the current-bias mode. 
The obtained values of the switching current, I = 
5.2 /l«A± 10 % ( J c ~ 10 /iA/ /mi 2 ) and the retrapping cur- 
rent, I r = 3 — 4.4 /iA, correspond to a small hysteresis of 
20 to 40 % and to an appreciable pair-current suppression 
factor 18], a = I a /I^ B ~ 0.8. To find I C AB w 6.4 pA we 
used an asymptotic junction resistance, i?N = 44 f2±10 % 
(corresponding to a value of p « 22 SI x /im 2 ). Us- 
ing a typical values of the specific junction capacitance, 
c w 50 — 75fF//tm 2 , we evaluate the McCumber pa- 
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FIG. 2: (a) Typical current voltage characteristics of a two- 
junction Al-SQUID, demonstrating a small remaining hys- 
teresis of the supercurrent. When the supercurrent is sup- 
pressed by magnetic field, BS — $o(?i + 1/2), where $o = 
2.07x 10 -15 Wb is the flux quantum, nonlinearity of the quasi- 
particle branch with 77 ~ 16 is observed, (b) TJM-model sim- 
ulation of an SIS-junction with the value of f3 c giving the best 
agreement with experiment. The model does not take into 
account the junction overheating effects, which were observed 
in experiment in form of a local back-bending of ZU-curves. 



rameter, (3 C ~ 0.8 — 1.2. In the junctions with higher 
transparency (with p = 10 — 12 x /tm 2 ), formed be- 
tween the first and the second shadows of Al, we found 
I = 10-13 /iA, J r = 7-9/iA, and 7 C AB « 13 /iA. Despite 
the smaller value of j3 c ~ 0.5 , a width of hysteresis was 
comparable to that of the lower-transparency junctions. 

A width of hysteresis of ~ 20 % found in most of the 
junctions agrees roughly with the RSJN model using the 
values of [3 C ~ 0.5 — 1 and an effective quasiparticle re- 
sistance at V < V c , Rj ~ (2 — 4)i?N (see, e.g., Fig. 4 
in Ref. 19]). However, more accurate calculations in the 
frame of the tunnel junction microscopic (TJM) model 
(see, e.g., Ref. and the references therein), accounting 
for the experimental value of a, produced the /U-curves, 
see Fig. EJb), in agreement with experiment, but corre- 
sponding to a lower value of f3 c w 0.3. The lower effective 
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value of f3 c indicates a stronger damping in our experi- 
ment than the predicted intrinsic damping due to quasi- 
particle tunnelling only. An additional damping mech- 
anism possibly arises due to a high-frequency shunting 
effect of the junction by the on-chip wiring impedance. 
This effect should be appreciable due to the relatively 
high values of i?N which are comparable to a typical 
microstrip-line impedance, R c ~ fO — 100 f2 (see, e.g., 
[13]), and can explain a similar hysteresis in the junctions 
of different transparency. An opposite effect of wiring 
was a larger hysteresis, up to 44 %, in those junctions 
which had the shortest connecting paths, ~ 100 /im long, 
to the wide leads of large capacitance. In practice, this 
contribution can be considerable in the systems of sub- 
micron tunnel junctions, but it can be reduced using, for 
instance, adjacent biasing resistors. 

The measurably larger retrapping current values can 
also be an effect of thermal fluctuations, activated pre- 
dominately in the resistive (R) state due to the junction 
self-heating. Using the model for noise-induced R^S 
transitions, Ref. [2l|, and the parameters of our junc- 
tions, we found a realistic estimate for the effective tem- 
perature in the R-state, T* ~ I K < T c . This estimate 
is consistent with the pronounced back-bending which is 
seen in the TV-curves in Fig. E^a) (cf., e.g., Ref. [12] )j 
and the observed ~ 20% reduction in the effective value 
of Vg at I > I Q . 

The quasiparticle branches of the JV-curves, like that 
shown by the dashed line in Fig. were measured in 
the SQUIDs with a supercurrent suppressed by a mag- 
netic field. The subgap conductance was found to be 
low, [G(0)] 1 ~ I kfl, which would allow low decoherence 
rates |4| in all-Al RSFQ-qubit circuits. High values were 
obtained for the nonlinearity parameter rj, which were 
larger than both in the Nb-based junctions 0, 0] and 
in the Al SIN-junctions 0. Similar to the case of SIN- 
j unctions, the values of r\ were found to depend on the 
junction configuration and size. In particular, the largest 
value of r\ « 32, obtained in the four-junction SQUID, 
exceeded ry w 16, measured in the two-junction SQUID 
with the two times larger junctions. In these two devices 
(not shown), the junctions were simply formed as over- 
laps of the straight perpendicular microstrips. A lower 
value of 77 = 8.8 was, however, found in a structure shown 
in Fig. with one of electrodes made of T" -shape (both 
right electrodes of the junctions; we call this configura- 
tion an " electron confinement shape" ) . The observed de- 
pendencies can be interpreted following the argument de- 
veloped in Ref. [^ for describing an enhancement of the 
Andreev reflection processes due to coherent two-electron 
(quasiparticle) diffusion in the junction electrodes. (See 
for a more detailed discussion Ref. Q.) 

To conclude, we investigated the damping in self- 
shunted Josephson junctions of type A1/A10 X /A1 of high 
transparency, fabricated using a shadow-evaporation 
technique. The /^-characteristics were featured by a 



critical supercurrent of the order of 10 \iA and a small 
hysteresis down to 20 %. Strong subgap nonlinearity of 
the quasiparticle branches was observed, which enables 
low decoherence in RSFQ-qubit integrated circuitry on 
the basis of Al junctions. 
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projects RSFQubit and EuroSQIP. 



* Electronic mail: Sergey.Lotkhov@ptb.de 
[1] K. K. Likharev, Dynamics of Josephson Junctions and 

Circuits (Gordon and Breach, New York, 1986). 
[2] K. K. Likharev, and V. K. Semenov, IEEE Trans. Appl. 

Supercond. 1, 3 (1991). 
[3] V. K. Semenov and D. V. Averin, IEEE Trans. Appl. 

Supercond. 13, 960 (2003); A. M. Savin, J. P. Pekola, 

D. V. Averin and V. K. Semenov, J. Appl. Phys. 99, 

084501 (2006b 

[4] F. Chiarello, cond-mat/0602464 

[5] A. B. Zorin, M. I. Khabipov, D. V. Balashov, R. Dolata, 

F. -I. Buchholz, and J. Niemeyer, Appl. Phys. Lett. 86, 
032501 (2005). 

[6] A. B. Zorin, E. M. Tolkach eva, M. I. Khabipo v, F.- 

I. Buchholz, and J. Niemeyer, cond-mat/0512615 
[7] S. V. Lotkhov, D. V. Balash ov7"Nl. 1- Khabip ov, F.- 

I. Buchholz, and A. B. Zorin, cond-mat/0605237 
[8] J. Hassel, H. Seppa, and P. Helisto, |cond-mat /0510189 
[9] A. J. Dahm, A. Devenstein, EX N! Langenberg, 
W. H. Parker, D. Rogovin, and D. J. Scalapino, Phys. 
Rev. Lett. 22, 1416 (1969). 
[10] G. J. Dolan, Appl. Phys. Lett. 31, 337 (1977); 

J. Niemeyer, PTB-Mitt. 84, 251 (1974). 
[11] E. L. Hu, R. E. Howard, L. D. Jackel, L. A. Fetter, and 
D. M. Tennant, IEEE Trans. Electron. Dev. 27, 2030 
(1980); L. D. Jackel, E. L. Hu, R. E. Howard, L. A. Fetter, 
and D. M. Tennant, IEEE Trans. Magn. 17, 295 (1981). 
[12] Yu. Makhlin, G. Schon, and A. Schnirmann, Rev. Mod. 

Phys. 73, 357 (2001). 
[13] V. Patel and J. E. Lukens, IEEE Trans. Appl. Supercond. 
9, 3247 (1999). 

[14] W. Chen, A. V. Rylyakov, V. Patel, J. E. Lukens, and 

K. K. Likharev, Appl. Phys. Lett. 73, 2817 (1998). 
[15] A. W. Kleinsasser, R. E. Miller, W. H. Mallison, and 

G. B. Arnold, Phys. Rev. Lett. 72, 1738 (1994). 

[16] R. E. Miller, W. H. Mallison, A. W. Kleinsasser, 
K. A. Delin, and E. M. Macedo, Appl. Phys. Lett. 63, 
1423 (1993). 

[17] Y. Naveh, V. Patel, D. V. Averin, K. K. Likharev, and 

J. E. Lukens, Phys. Rev. Lett. 85, 5404 (2000). 
[18] A. B. Zorin, K. K. Likharev and S. I. Turovets, IEEE 

Trans. Magn. 19, 629 (1983). 
[19] D. E. Prober, S. E. G. Slusky, R. W. Henry, L. D. Jackel, 

J. Appl. Phys. 52, 4145 (1981). 
[20] T. Hoist, D. Esteve, C. Urbina, M. H. Devoret, Physica 

B, 203, 397 (1994). 
[21] E. Ben- Jacob, D. J. Bergman, B. J. Matkowsky, 

andZ. Schuss, Phys. Rev. A26, 2805 (1982). 
[22] P. N. Dmitriev, et al., IEEE Trans. Appl. Supercond. 13, 

107 (2003). 

[23] F. W. J. Hekking and Yu. V. Nazarov, Phys. Rev. Lett. 



4 



71, 1625 (1993); F. W. J. Hekking and Yu. V. Nazarov, 
Phys. Rev. B. 49, 6847 (1994). 



